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Abstract 

We consider generalized Forchheimer flows of either isentropic gases or slightly compress¬ 
ible fluids in porous media. By using Muskat’s and Ward’s general form of the Forchheimer 
equations, we describe the fluid dynamics by a doubly nonlinear parabolic equation for the ap¬ 
propriately defined pseudo-pressure. The volumetric flux boundary condition is converted to a 
time-dependent Robin-type boundary condition for this pseudo-pressure. We study the corre¬ 
sponding initial boundary value problem, and estimate the and (with 0 < o < 1) 

norms for the solution on the entire domain in terms of the initial and boundary data. It 
is carried out by using a suitable trace theorem and an appropriate modification of Moser’s 
iteration. 


1 Introduction 

The most common equation to describe fluid flows in porous media is the Darcy law 

-Vp=^v, (1.1) 

where p, v, p, k are, respectively (resp.), the pressure, velocity, absolute viscosity and permeability. 

However, this linear equation is not valid in many situations, particularly, when the Reynolds 
number increases, see wm- Even in the early work, Darcy [Ij already acknowledged the deviations 
from equation There have been many investigations into what equations for hydrodynamics in 

porous media to replace Darcy’s law O), see [3l l25l[26ll28ll33j and references therein. Forchheimer 
OE] established the following three nonlinear empirical models: two-term Forchheimer equation 

— Vp = av + b\v\v, (1-2) 


three-term Forchheimer equation 

— Vp = av + b\v\v + c\v\‘^v, 


(1.3) 


and Forchheimer’s power law 

— Vp = av + d\v\‘^~’'v, for some real number m G (1, 2). 


(1.4) 
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Above, the positive constants a, b, c, d are obtained from experiments. 

While mathematics of Darcy’s flows have been studied intensively for a long time with vast 
literature, see e.g. m, there is a much smaller number of mathematical papers on Forchheimer 
flows and they appeared much later. Among those, there are even fewer papers dedicated to 
compressible fluids. (See |28j and references there in.) 

In order to cover general nonlinear flows in porous media formulated from experiments, gener¬ 
alized Forchheimer equations were proposed. They extend the models (jl.2l) - (ll.4p and are of the 
form 

N 

— Vp = ^ (1.5) 

i=0 

These equations are analyzed numerically in [7lll9p 27j , theoretically in [2 l[T0l - [T^ [T6] for single-phase 
flows, and also in |14lll5j for two-phase flows. 

Our previous analysis piinHiaiis] was focused on a simplified model for slightly compressible 
fluids. Though such a minor simplification is commonly used in reservoir engineering, the mathe¬ 
matical rigor is compromised. Furthermore, since the model does not specify the dependence on the 
density, its applications to gaseous flows would be inaccurate and might present artificial technical 
difficulties. The goals of this paper are: (a) Developing a more accurate model for generalized 
Forccheimer equations for gases, and (b) Analyzing it without making any simplifications. 

For goal (a), we first have to modify (II.5p to reflect the dependence on the density. We return 
to an idea by Muskat and Ward. By using dimension analysis, Muskat [25] and then Ward [33] 
proposed the following equation for both laminar and turbulent flows in porous media: 

— Vp = p °‘~^where / is a function of one variable. (1.6) 


In particular, when 0 = 1,2, Ward |33| established from experimental data that 

— Vp = -|- cf-^\v\v, where cp > 0. 

k Vfc 


(1.7) 


Combining (jl.5ll with the suggestive form (|1.6p for the dependence on p and v, we propose the 
following equation 


N 


— Vp = ^ (1.8) 

i=0 

where A" > 1, oq = 0 < Q^i < ... < ow are real numbers, the coefficients oq, ..., aw are positive. 

Here, the viscosity and permeability are considered constant and we do not specify the depen¬ 
dence of Oj’s on them. Our mathematical exposition below will allow all Oj > 1 in (11.81) . In practice, 
we can simply take ow < 2 in (|1.8p or use the popular model (jl.7p . Even in these cases, the results 
obtained in this paper are still new. 

Multiplying both sides of (|1.8h by p gives 


g{\pv\)pv = -pVp, 

where g : M"*" —>■ M"*" is a generalized polynomial with positive coefficients defined by 


(1.9) 


g{s) = aos°‘° -k ais"^ H--k aws”^ for s > 0. (1-10) 

We will study the following two types of compressible fluids. 

1. Isentropic gases. For isentropic gases, the constitutive law is 


p = cp'^ for some c, 7 > 0. 


(1.11) 
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Then from (11.91) and (jl.llll follows 

g{\pv\)pv = —pS/p = —Vu with u =-. 

7+1 

Solving for pv from this equation yields 

pv = —K{\Vu\)Vu, 

where the function K : M'*' —)• M'*' is defined for ^ > 0 by 

K{^) = ^ with s = s(^) being the unique non-negative solution of sg{s) = 


Recall the continuity equation 

(j)pt + div{pv) = 0, 

where constant (p £ (0,1) is the porosity. Rewrite 

,7 + 1 , ^ 

' ' \ ^,_L1 


1 


p=( -It+iu'' with a =--€( 0 , 1 ). 

^ ^ c-y ^ 7 + 1 


C7 ' 7 + 

Combining (jl.l5p and (jl.l3|) with relation (I1.16|) . we have 

1 / C7 


= V-(R:(|Vn|)Vu)). 


2. Slightly compressible fluids. The equation of state is 

1 dp 1 

= const. > 0 . 

p dp K 

Note from (11.181) that pVp = kV p. Then combining this with (jl.Op gives 


( 1 . 12 ) 

(1.13) 

(1.14) 

(1.15) 

(1.16) 

(1.17) 

(1.18) 


g{\pv\)pv = —Vu with u = np. (1-19) 

which is the same as (I1.12p . Thus, we obtain formula ()1.13p for pv. By combining (I1.13P and (11.151) 
we have 

ut='^V ■{Ki\Vu\)Vu). (1.20) 

(p 

Observe that we can write ut = {u^)t with A = 1 in (jl.20l) . Therefore, the two equations (I1.17P 

and pi.20p are the same except for the factors on the right-hand sides. Then by scaling the time 

variable in both ()1.17p and (jl.20l) . we obtain, for both isentropic gases and slightly compressible 
fluids, the unified equation 

(n^)i = V-(R:(|Vn|)Vn)) with A G (0,1]. ( 1 . 21 ) 


This will be the partial differential equation (PDE) of our interest. It is derived and will be analyzed 
without any simplifications (goal (b) above). Although the case of isentropic gases, i.e. A < 1, is 
the main focus, the analysis will also cover the case of slightly compressible fluids, i.e. A = 1, at no 
extra cost. 

In case of ideal gases, i.e., 7 = 1, we can derive from (|1.2ip a PDE for pseudo-pressure 

7 + 1 

In general, we rewrite u in (|1.12p as u = dp t" for some d > 0, hence it is approximately the 
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pseudo-pressure for isentropic gases (11.111) . For simplicity, we refer to u as the pseudo-pressure. 
Therefore, equation (11.211) is a PDF describing the dynamics of the pseudo-pressure u. 

Regarding the boundary conditions, we focus on the volumetric flux condition v ■ V = if), which 
results in pv ■ V = ijjp, or, 

— K{\Vu\)Vu ■ P = ipu^, (1-22) 

where A G (0,1), </? = ™ case isentropic gases, and \ = 1, g) = 'll )/k in case of slightly 

C'y 

compressible fluids. Here, P is the outward normal vector on the boundary. 

From mathematical point of view, equation (jl.2ip for A < 1 is a doubly nonlinear parabolic 
equation, which is an interesting topic of its own. Research on doubly nonlinear parabolic equations 
follows the development of general parabolic equations [2Tl[22j and degenerate/singular parabolic 
equations 016 ]. However, it requires much more complicated techniques. See monograph 
review paper [18] and references therein. For other developments, see e.g. [Tl l^l^l30ll32] . 

There are two issues that did not attract attention in most existing papers: (1) Robin boundary 
condition, and (2) Estimates for super-critical case. For instance, [2^1^ give global L°°-estimates 
but for homogeneous Dirichlet boundary condition in the sub-critical case, see discussion in Remark 
[Q below. Regarding interior estimates, the common result is 

ll'“llL°°(Q{iJ/ 2 ,r/ 2 )) < Li^(Q(R^T)) ■*" ^)’ (1-23) 

where q G (l,oo), Cr^t > 0 depends explicitly on R and T, while Q{R,T) denotes the cylinder 
Br X (—T, 0). Surnachev [29] improves it to 

lkllL-(Q(R/2,T/2)) < a{q),l3{q) > 0. (1.24) 

We call (|1.24p a quasi-homogeneous estimate (with respect to ||if||L 9 (Q(_R,T))). The global version 
of (I1.24P is not known for doubly nonlinear equations, though it was established for degenerate 
equations, see e.g. m- 

In this paper, we focus on both topics (1) and (2) listed above. In this case, the boundary 
condition (ll. 22 p gives rise to high power boundary integral which cannot be treated by the standard 
trace theorem. Therefore, we derive and utilize a new, suitable trace inequality to obtain bounds 
for the solutions of (|1.21l) in terms of initial and boundary data. For L°°-estimates, we make 
some technical improvements in order to overcome the non-homogeneity of function K{-) and non¬ 
zero boundary data. We carefully modify Moser’s iteration |24j and obtain quasi-homogeneous 
estimates. Our results are for both (spatially) interior and global estimates, hence, extend the 
previous interior improvement (|1.24p . 

Throughout this paper, U is an open, bounded subset of M"", with n = 2,3,..., and has C^- 
boundary F = dU. For physics problems n = 2,3, but we consider here any natural number n > 2. 
Hereafter, we fix the functions g{s) in (11.91) and (jl.lOl) . Therefore, the exponents a* and coefficients 
ai are all fixed, and so is the function K{^) in (11.141) . Also, our calculations frequently use the 
following exponent 

a = ^^^(0,1). (1.25) 

We consider the initial boundary value problem associated with ()1.21l) and (|1.22p . specifically, 

= V-(A:(|Vu|)Vu) inC/x(0,oo), 

< tt(x,0) = uq{x) in U, 

^K{\Vu\)^ + ipu^ = 0 onFx( 0 ,oo), 


(1.26) 
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where uo(x) and are given initial and boundary data, respectively. Again, F denotes the 

outward normal vector on F. Here, A is a fixed number in (0,1] for the remaining of the paper. 

The current article is focused on studying non-negative solutions of problem (|1.26p . Section 
[2] contains new trace theorems and multiplicative Sobolev’s inequalities, which are suitable to the 
Robin-type boundary condition (ll.22p . as well as the nature of our equation’s double nonlinearity. 
In section O we estimate L“-norms of the solutions for all a > 0, in terms of initial and boundary 
data. These will also be used for later gradient and L°° estimates. In section 01 we present estimates 
for the gradient’s L^““-norm for time t > 0. In section (H we estimate the L°°-norm of the solution 
in any compact subsets of the domain. Due to the basic Lebesgue norm relation in Proposition 15.21 
the Moser’s iteration is of a non-homogeneous form (I5.22p . We deal with this by using Lemma lA.21 
and obtain in Theorem 15.,11 the quasi-homogeneous estimate. Section [ 6 ] is focused on estimating 
the solution’s L°°-norm on the entire domain. In Proposition 16.21 the L'^“-norm is bounded by 
the L"+^i -norm, where, unlike the interior case, both k > 1 and > 0 depend on a, and the 
constant depends on the boundary data. To manage the powers during iterations, we construct in 
Lemma [6~4l two controlling sequences (aj}'^Q and (/3j)^o- Using these sequences for iterations, we 
obtain in Theorem 16.61 the quasi-homogeneous estimates for the L°°-norm, and in Theorem 16.81 the 
ultimate estimates in terms of initial and boundary data. The Appendix contains key Lemma lA.21 
in implementing Moser’s iteration with non-homogeneous inequalities. 


2 Auxiliaries 


First, we recall elementary inequalities that will be used frequently. Let x,y > 0. Then 


(x + y)^ < 2P(x^ + yP) for all p > 0 , ( 2 . 1 ) 

(x yY < x^ + y^ for all 0 < p < 1 , ( 2 . 2 ) 

{x + yY <{x^ + y^) for all p > 1, (2.3) 

x^ < x°‘ + x'^ for all 0 < a < /3 < 7 , (2.4) 

particularly, 

x^ <1 + x"^ for all 0 < /3 < 7 . (2.5) 


Second, we establish particular Poincare-Sobolev inequality and trace theorem for studying our 
doubly nonlinear equation with Robin-type boundary condition. 

For any 1 < p < n, we denote by p* its Sobolev conjugate exponent, that is, p* = Y^- 

Lemma 2.1. In the following statements, u(x) is a function defined on U. 

(i) If ot> s>f), a>l, and p > 1, then for any |u|“ E W^’^{U) and e > 0 one has 

f \uYda<£ f \uY~^\VuYdx + Cl f \uYdx + {c2a)p^e~p^ f \uY~^Fodx, ( 2 . 6 ) 

Jr Ju Ju Ju 

where ci,C 2 > 0 are constants depending on U, hut not on u{x),a, s,p. 

(ii) lfn>p>l,r>0, a>s>0 , a> and a > !l(!±£zlii^ then for any e > 0 one has 


'U 


f 9 9{ot — s+'p) 

|"+'’dx<e / \uY-^\S/uYdx + £~^2 {c^m) 

Ju 




+ 2®(“-"+P)c4^|t/| 


6>(ct(p—l)+s—p) 


\u 


\a+r 

IL“ ’ 


(2.7) 
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for all E W^’P{U), where 

a — s+p rn r + 9(s — p) 

m = -, 0 = — -—-, pi = -^--- 

p n[p — s) + ap \ — 0 

and constants 03,04 > 0 depend on U,p, but not on u{x),a, s. 

(iii) If n > p > 1, a > s > p, and a > , then for any e > 0, one has 


( 2 . 8 ) 


f \u\°‘da< 2 e f \u\°‘ ^\Vu\^dx + ci\\u\\fa 

Jr Ju 

+ e-^L»i,„||n||2r (2.9) 

for all functions u{x) satisfying \u\°‘ E W^'^{U) and \u\"^ E W^’'^{U), where m is defined by 

S ( 2 - 10 ) 




= 2®(“-^+?’)(o2a)i^of |?7| 


0(a(p—l)+s—p) _ 6{a — s+p) p 1 _ 0p _ 

^ ) -^2,0 = 2 1-® (o2a)p-i 1-® (03m) i-«. (2.11) 


Proof. We make a couple of comments before starting the proof. First, the boundary integrals in 
(12.6p and ()2.9p are in the sense of traces of |u|“ on F. Second, observe that the conditions on u{x) 
do not guarantee that the right-hand sides of p2.6h . (j2.7p and (j2.9p are finite. In case they are not, 
these inequalities are understood to be trivially true. Therefore, the following proof only needs to 
cover the case when those right-hand sides are hnite. 

(i) We recall the trace theorem 

/ < Cl / \(f\dx + C 2 / \V(f)\dx, 

Jr Ju Ju 

for all (f E W^’^{U), where ci and C2 are positive constants depending on U. Applying this trace 
theorem to ()> = |n|", we have 


u\^da<ci / |u|“dx-|-C2Q! / |u|“ ^\Vu\dx 

Ju Ju 


( 2 . 12 ) 


/ Q —s , (p—IJQ+S—p 

Rewriting C2a|u|"“^|Vu| in the last integral as a product of p |Vu| and C 20 Le~^''^u p , 

and applying Young’s inequality with exponent p and p/{p — 1), we obtain inequality (|2.6p . 

(ii) Since a > s, the number m defined by (|2.8I) is greater or equal to 1. Then applying 
Sobolev-Poincare inequality to Ini™ yields 


Ilp* < C 3 ||V(|n|™)||i,PC4 / |n|™dx 

Ju 


(2.13) 


where C3 and C4 are positive constants depending on U and p. Note that by definition 
we have (m — l)p = a — s. Hence (I2.13p can be written as 


of m. 


i/p’ 

|nr™dx) < 


u 


C 3 m(^ j |n|" *|Vn|^dx^ ^ + C 4 J \u\'^dx. 
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Raising both sides to the power 1/m < 1 and using inequality (|2.2I) . we obtain 




u^dx 


Note that 


q = p m = 


n{a — s + p) 


n — p 


Then (|2.14|) yields 


\u\\l<i < 




+ cT \\u\\l^ 


(2.14) 

(2.15) 

(2.16) 


Since n > p, r > 0, and a > one has a < a + r < q. Then 


1 


where 0o € (0,1) is defined by 


0n = 


a + r q 


rq 


_ 00 , 1 - 00 

I ? 


a 


{a + r){q - a) 
Then interpolation inequality and (|2.16l) give 


with q = p*m, 


(2.17) 






I 


+ cT\\u\\l- 




Raising both sides to the power a + r and applying inequality (12.Ih for exponent 6 o{a + r) yield 

+ ea{a + r) 

71 


u 


lu 


a —S + p 


dx < 2®“^“^'’^|(c3m) ^ J It" 

= J u^-^\Vu\Pdx 

I 9 ^ 0 II ||(l-6»o)(a+r) n ||0o(a+r) 


\u\ 


6»0(«+?•) 1 II II (l-0o)(Q:+r) 


u 


where 


0 = Mei±h which is the same as in 


(2.18) 

___ (2.19) 

a — s + p 

Since a > "^F+^-v) , we have 0 E (0,1). Then applying Young’s inequality to the first term on the 
right-hand side of (I2.18P with powers ^ and 7^, we obtain 

[ I /* 0 gn(a + r) gn(Q; + r) j (l — Qo)((^-F^)—L— 

/ u°‘^^dx<e / tt" ^\Vu\^dx + s ^-<^2 i-e (csm) i_e 

Ju Ju 


_l_ 2®o(a+?-)o 


fln(a; + r) 


m 




Since a > ^7; then m < a. By applying Holder’s inequality to bound the L^-norm of u on 


the right-hand side of (I2.20p by |f7| m a ||M||i,a, we obtain 


u\°'^'^dx<e I |u|" ^{Vul^dx + e i-e 2 i-» (cam) 


9 6>n(a; + r) 


'C/ 


epCa+r) 1 (l_0p)(Q,_|_r)-7 


\u\ 


lu 


t>o(a; + r) 


+ 2 ®o("+^)c 4 (2.21) 
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Re-calculations of the powers: 

0o(a + r) = 6(a — s + p), 6Q{a + r)/m = Op, 




1 — 6 \ a + r J I — 6 

r + 9{s — p) 

= a-\ - 1 — 0 - = o: + Pi, 


1-0 


eoia + r){— - -) = 0{a - s + p)( 
m a \ 


P 


1 \ 0{a{p — 1) -|- s — p) 


a — s + p a 


a 


Thus, inequality ()2.7I1 follows ()2.21l) . 
(iii) Define 


r = 


s — p 


p — 1 


( 2 . 22 ) 


Given e > 0. First, we apply inequality (12.61) . and then estimate the last integral Jjj \u\^~^'^dx 
in (j2.6|) by using (12.7p with the parameter s in (j2.7p being set as 


p 1 


e{c2a) p-iep-i =ep-i(c2Q;) p-i = (e ^C2a) p-R 


This results in 


[ u"d(j<2e [ u'^-^\Vu\Pdx + ci [ u^dx 

Jv Ju Ju 


p 1 p 6 9(a+s—p) 9p rvJ-i 

+ {c2a)p-^e p-i (e“^C2a)p-i i-® 2 ^ (csm) 


+ e“ 7^ (caa) 7^ 2®(“-"+P)cf 117| 


9{a{p—l) + s—p) 


I a+r 

IL" ’ 


where 0 is defined in (j2.8p . which is the same as in ()2.10l) since r is now specified by (j2.22l) . 
Therefore (12.9p follows. Finally, we check the conditions on the exponents in (i) and (ii) to validate 
our calculations. Since a > s > p > 1, we have r > 0, {p — s)/{p — 1) < 0, and only need to check 
a > !i(l±£zP)^ With r defined by (j2.22p . this, in fact, is a > n{s — p)l{p — 1), which is already one 
of the assumptions on a. The proof is complete. □ 

In our particular case, we have following lemma. Define 


(5 = 1 — A G [0,1), a* = n(a — (5)/(2 — a) and pQ =- 


a — 6 


— a 


(2.23) 


Lemma 2.2. Assume a>S,a>2 — 6 and a > npQ. Let c* = max{ci, 02 , 03 , 04 } with oi, 02 , 03,04 
in Lemm,a\2.1\. and 


0 = 0 ^ 


1 


(1 — a)(a/a* — 1) 


G(0,1). 


(2.24) 


Then one has for any e > 0 that 

J |u|“d(T < 2e J \u\°‘~^^~‘^\Vu\‘^~‘^dx + c^,\\u\\j^a(^ij^ 


4 - On I TJ ||„.||«+/^1 

' \\^\\La(U 


(2.25) 
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where 


def /io(l + ^(1 - a)) def 

hi — hi,a — -~ ~ 

(2-a)(l + e(l-a)) 

Dg „ = a^\U\ 

6(^OL+5—a) (2 —a)(l+0(l —a)) 

= 2 i-e (c*a) . 


+ 


0(2-a) 


1 — a (1 — 0)(1 — a) ’ 


(l-a)(Q + /Jn)9 


(2.26) 

(2.27) 

(2.28) 


Proof. We apply inequality (12.9p in Lemma[2)T]to p = 2 — a and s = 2 —5. We recalculate exponents 
for these particular values. Note, a — s + p = a + 6 — a. Then from (12.8p and (|2.22p . 


m = 


a + (i — a 
2 — a 


r = 


s — p 

p — 1 


a — 6 
1 — a 


= ho- 


Also, 9 in (|2.10p becomes (I2.24p . the number pi in (j2.8l) is the same as in (I2.26p . The exponent 
of \U\ is 0(a(l — a) + a — 6)/a = 9{a{l — a) + (1 — a)po)/a = (1 — a)6{a + po)/a. The exponent 
of e~^ in the last term of (12.9p is p 2 - Also using the fact ci, 02,03,04 < 0* and m < a, we have 
Di,a < 1^3,0 and Z?2,a < D 4 ,a- Therefore, we obtain (I2.25P from ()2.9I1 . □ 


Next is a parabolic multiplicative Sobolev inequality. 


Lemma 2.3. Assume 


a >2 — 5 and a > a*. 


(2.29) 


// T > 0, then 
rT 

lO JU 


L 


< ( 050 ^"“)-( / / \uY+^-^dxdt+ I I \uY+^-^\Vu\‘^-'^dxdt 


o+5-2|^„,|2-a 7 i.\ “+'5- 


0 JU 


10 JU 


sup ( / \u{x,t)Ydx 
te[o,T] ^ Ju 


1-fl 

\ “ 


where 05 > 1 is independent of a and T, and 


9 = 9 a = 


n{a+S—a) 

In ease U = Br - a ball of radius R - one has 
T 


L 


0 J Bu 




fl 

'-Jo J B 


uY+^-^^dxdtP [ [ \uY+^-^\Vu\‘^-‘^dxdt 

Jo J Bf 


a + 5 — a 


(2.30) 


1- / \ def .. 2 a a 6 -it/ /'oQif 

- , , , n = K{a) = 1 +-= 1 + (a - 6 ){ -). (2.31) 

X _j_ a{^-a) n a Of* a 

' ft ( /-v_L A_/I 


sup ( / \u{x,t)\'^dx] , (2.32) 

te[o,r] ^Jbr 


where cg > 1 independent of a, R and T. 
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Proof. Note that by definition of k and 9 we have ^ ^ where po = ■ Then 

interpolation inequality gives 




u 


KO. 


\u\^°dx] ■ 


u 


lul^dx 


u 


(2.33) 


We estimate the first integral on the right-hand side of (j2.33l) . We recall the following Sobolev’s 
inequality 


2-a 


u;||^(2-a)* <C7( / \wg °'dx+ / iVtcf °‘dx 


u 


lu 


where (2 — a)* is the Sobolev conjugate of 2 — a, and cy > 1 is independent of a. Applying this 
inequality to re = lul™ with m = > 1, we obtain 


n(Q; + 5 —a) 

|u| dx 


u 


(2-a)* 


< C7 


U 




l"+<5-2|w.,,|2-a I \ 2 -a 


U 


lU 




(2.34) 


Note that Ka9 = a + 6 — a and (2 — a)*/po = {2 — a)/{a + 6 — a). Then raising both sides of 
inequality (12.341) by (2 — a)*Ka6/po = 2 — a results in 


u 


\urdx] < 4-“ 


\u 


-b m 


u 


f \u 

Ju 


a+5-2|Y7,,|2-a 


\Vu\^-^dx 


(2.35) 


{1—9)k 


Raising both sides of (I2.33P to the power na, and using inequality (I2.35p . we get 

J \u\^^dx <{c7mf-^(^j |u|“+^-“dx + m2-“y' \u\'^+^-^\Vu\‘^-‘^dx'^ J |u|"dx) 

Integrating this inequality in t from 0 to T and taking the supremum of the last integral for 
t E [0, T] yield 


r/ 

/o Ju 


uy^^'dxdt < {cTrnf '^(^J J 

sup ( / |u(x,t)l"dx 
te[o,T] ^ Ju 


a-\-S—a 


dxdt + 


fl 


10+5—21 v7oi 12—a 


|Vuf °‘dxdt 


{i-e)K 


Taking both sides to the power ^ , we have 


0 JU 


KOL 


u\'^^dxdt] < (crm) 


2 —a 

KOt 


L 


u 




sup 
te[o,T] ''Ju 


0 Ju 

u{x, t)\^dx 


fl 


l“+'5-2|w,,|2-a 


\Vuy-^dxdt 


a-l-i5 —a 


Note that m < a, then we obtain ()2.30p with C5 = Cy “. 

Consider the case U = Br. The Sobolev inequality for Br is 




-1 


w\^-^dxY~’^+cA I \Vw\^-^dx 


Br. 


|2 — 2-a 


Br 


(2.36) 


where cs > 1 is independent of R. Repeating the above proof with C7 replaced by 03(1 + i? ^), we 
obtain (j2.30p with C5 replaced by [03(1 -|-R~^)]^~“. Therefore, we obtain (j2.32p with cq = Cg~“. □ 

It is noteworthy that the explicit exponents and constants in Lemmas 12.21 and 12.31 play an 
important role in Moser’s iterations in sections [5] and [6l 
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3 L“-estimates 


We start studying problem (jl.26jl . Hereafter, u(x,t) denotes a non-negative solution of (jl.26jl . 

Regarding the nonlinearity of the PDE in (ll.26p . we recall that the function has the 

following properties: it is decreasing in maps [0,oo) onto (0, ^], and 


di 

(i + 0 “ 


< < 


d2 

(1 + 0 “ 


d3(0"“ -1) < < ^20"“, 


(3.1) 

(3.2) 


where di,d 2 ,d 3 are positive constants depending on g(s). (See [2] for the proof.) 

Notation. The symbol C denotes a generic positive constant with varying values in different 
places, while Ci,C 2 ,--- and ci, C 2 ,... have their values fixed. The constants C, C, Cj, Cj, for i,j = 
1, 2, 3,... , used in calculations can, otherwise stated, implicitly depend on number A, function 
the space dimension n, and U. We also use the positive and negative parts notation /+ = max{/, 0} 
and f~ = max{—/,0}, and brief notation ||</5(t)||^oo(p) for ||</5(-,t)||x,cx)(p). 

Recall 5, a* and go are defined in (|2.23l) . We assume hereafter that 


a > 5, i.e., the number go is positive. (3.3) 

The case a < (5 is much simpler, see Remark 13.21 below. We begin with a differential inequality 
for fjj u°‘dx. 

Lemma 3.1. Assume a > 2 — <5 and a > ngo- Then 


^ J u‘‘(x,t)dx+Ci J^\Vu(x,t)\'‘ “u"+* '^(x,t)dx <0\\u{t)\\ltii,f+C\\<f> (i)|li~(r)l|t‘(t)IIS.((;) 

+ C||^V*)lllS,r)IIMOII^i('S, +C||»>-{«)|lI5|pi|a(i)|lS('‘t;), (3.4) 

with gi and 9 in Lemma l+il where the positive eonstants Ci = Ci^a and C > 0 depend on a. 
Consequently, there is C 2 = +' 2 , 0 - > 0 such that 

^ ^ u^dx + Cij^ I < C 2 (1 + Ilv^"(t)|| (1 + ^ u^dx) (3.5) 

Proof. Multiplying both sides of the first equation in (jl.26l) by , integrating over domain U, 

and using integration by parts, we have 


a dt 



V ■ {K{\Vu\)Vu)u‘^+^-^dx 

j 

(a-A) [ K{\Vu\)\Vufu^+^-^dx+ [ K{\Vu\)^u^+^-^da 
Ju Jr du 

(a-A) [ K{\Vu\)\Vu\‘^u'^+^-^dx- [ u^ipda. 

Ju Jr 


Using relation — 1) in (|3.2p . one has 


(3.6) 


- / K{\\/u\)\S/u\^u^+^-'^dx <-do / |Vrt|2-“n"+'5-2dx + d3 / u^^+^-^dx 


lu 


lu 


lu 


<-da / \Vu\^-^u^+^-^dx + C\\u\\lt 
Ju 


a^5—2 


(3.7) 
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Estimate the last integral of (j3.6p by using the trace theorem in Lemma 12.21 we have 


J u°^ipda < J 


— I u°‘ipda < J u'^if da < \\(p J u'^da 

(r){2e f |nr+'5-2|Vu|2-“dx + C||n||2«(c7) + 

J u 


< 


\L°°( 


(3.8) 


Combining (|3.6I) . (|3.7I) and (13.81) . we have 

u°‘dx <—{a — X)d3 f \Vu\‘^~°'u‘^~^^~‘^dx + C\\u\\1t^~‘^ + \\‘f~\\L°°{r) 


a dt 


lu 




lu 

a+(5-2|Y7„,|2-a 


|Vnp-“dx + C||u||2.(t;) + Ce-T^||u||“+g^) (3.9) 


HU)j- 


The case || 
Select £ = 


L°°(r) 

d?, 

4||¥’“llLo°(r) 


= 0, inequality (13.41) immediately follows (13.91) . Consider 
in (|3.9p . Then we obtain 


L°°(r) / 0 - 


<-(a- A-^)d3 ^|Vu|^ dx + CWuW’ltljj) + C\\(p ||Loc(r) ||u||2c.(f/) 

+ + (3-10) 


Note that ^2 + 1 = (^i_a)(i-e) ■ Multiplying both sides of (13.101) by a/X, we get inequality (|3.4I) . 

Now, we prove (13.5|) . Denote /3 = 1 + fJ-i/a. Note that, on the right-hand side of (13.41) . the 
maximum power of ||rt||L“ is = af3, and the maximum power of ||(^~||Loo(r) is By 

applying inequality (12.51) . 


I ||a+(5-2 II iia n no+Zio <1 , |L,||“/3 


\W llL°°(r), 

Therefore, it follows from (13.4p that 


’"llLV)<i + 


2-a 

-|| (l-a)(l-e) 

llL°°(r) 


A 

dt 


u^dx + Cl 


'u 


[ |Vn| 

Ju 


2—a^a+5—2 


dx < (7(1 -b 


’)(i + NIl^) 


iiL°°(r) 

<C(l + ||^-||'^^)(l + ||n||^)^. 

Thus, inequality (|3.5D follows. □ 

Remark 3.2. In case a <6, by combining US. 6\) and the simple trace inequality 112.6\) . we can find 
X d 


a dt 


[ u°dx= [ V ■ {K{\Vu\)Vu)u'^+^-^dx 

Ju Ju 

= -Ci J K{\Vu\)\Vu\‘^u^+^-^dx + C^(^j |n|"dx 


+ 


(3.11) 


where C^p> Q depends on the function ip. Since po < 0 in this case, we can apply Holder’s inequality 
and easily derive a differential inequality for JjjU°’dx, and consequently obtain its estimates. Hence, 
there is no need to use more involved trace inequality 112.9\} . which is essential in the proof of Lemma 
3.1\ for the case a > 6. Therefore, we refer to a > 6 as the super-critical case, and a < 6 as the 
sub-critical case. Note that papers fall into the latter case. 
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We now have local (in time) estimates for solutions. 
Theorem 3.3. Assume a > 2 — 6 and a > npo. 

(i) // T > 0 satisfies 


f\i + < ^(l + , (3.12) 

Hpf 

where = Cs^q, = C2[1\Iol > 0, then for all t E (0,T]; 

J^u°‘{x,t)dx < ^(l +J^Uo{x)dx^ -CsJ (1 + (3.13) 


(ii) Consequently, ifT>0 satisfies 


(1 + 


2 — a 

(l_a)(l_e) 

L°°(r) 


)dt < 


I _ 


1 + / UQ{x)dx 


'u 


-uih 


then 


<0 JU 

where C > 0 depends on a. 


J u°‘{x,t)dx < 2 ^ + j no(x)dx^ for all t ^ 

J J u°'~^^~‘^{x,t)\Vu{x,t)\‘^~°'dxdt < c(^l + J UQ{x)dx'^, 


Proof, (i) Let ^(t) = 1 + j^u‘^{x,t)dx. Then inequality (|3.5I) can be written as 

jV{t) < C2(l + 


Solving this differential inequality gives 


V{t) < {y(0)-^i/“ - ^ j\l + \\<p-{r)\\^^^)dT] 

for all t € (0,T], with T > 0 satisfying (|3.12p . Hence we have (|3.13p . 

(ii) When T > 0 satisfies ()3.14p . inequality ()3.15l) easily follows ()3.13l) . 
Integrating inequality ()3.5p in time and using (I3.15p . we have 


(3.14) 


(3.15) 

(3.16) 


(3.17) 


[ [ u"+'5-2(x,t)|Vu(x,t)p-“dxdt 
Jo Ju 

< C y no(x)dx + C' J (^1 + J u°'{x,f)dx^ 


i+ai/o- 


1 + 


(l_a)(l-e) 1 ,, 

L°°(r) 


< 


C J UQ{x)dx + c(^l + J UQ{x)dx^ ^ J 


1 + 


2-a 

(l_a)(l-0) 

L“(r) 


dt. 


Combining this with the bound ()3.14p for the last integral, we obtain (|3.16p . 


□ 
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4 Gradient estimates 


In this section, we estimate Jjj |Vttp °'{x,t)dx for t > 0. Same as in [2], we will use the following 
function 

H{^) = / K{^)ds for ^ > 0. (4.1) 

Jo 

The function H{^) can be compared with ^ and K(^) by 

KiOf < H{^) < 2Km\ 

and hence, as a consequence of (|3.2p and ()4.2p . we have 

- 1) < H{i) < 2d2l'-“. 

Proposition 4.1. Assume 


a > max{n/io, A + 1 + /io}- 


(4.2) 

(4.3) 

(4.4) 


If t > 0 then 

f! 


^\{u^)t\‘^dxdT + f |Vrt(x,t)p “dx 

Ju 

+ CJC{t) + C (l + ||vp"(t)||J^^^^)(1 + ||u(T)||"+J)))dT, (4.5) 


<CTo 

where C > 0 depends on a, 


£o = I UQ{x)dx+ I |VMo(a^)|^ “dx + / UQ’''^(x)<y9’''(x, 0)d(T, 


'u 


lU 


JC(t) — 1 + 


M3 


+ 


with 


def 

h3 = h3,a = 


n |v?t(x,r)|“-^-id(Tdr, 

' 

(2 — a)a 


(4.6) 

(4.7) 


(1 — a)(a — (A + 1 + Ho)) 


Proof. Multiplying both sides of the PDE in (I1.26P by Ut = integrating over U and 

using the boundary condition, we obtain 

1 d 


1 

A 




'u 


dd(|Vu(x, t)|)dx = — y ipu^utda 


1 d f 

1 dt dp 


A + 


u^~^^(fda + 


u^^^(ftdcr. 


(4.8) 


A + 1 


Multiplying the equation by A + 1, and applying Young’s inequality to the last integral yield 


A + 1 


»■-"!+),p<ix + 4( 


d / A + 1 


lu 


dd(|Vu(x,t)|)dx + J u^~^^(pdaj = J u^^^iptda 

< J u°^da + J a-A-i dcr, (4.9) 
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To estimate the second to last boundary integral, we use the trace inequality (|2.25p in Lemma l2.2l 

f \u\‘^da<2e f + 

t/ r 'z u 


for any e > 0. Using this inequality in (I4.9p . we have 


A + 1 


u^-\u^)ldx + ^{ 


lu 


d / X + 1 


dt\ 2 

< 2s \u 

Ju 


'U 


H{\Vu{x,t)\)dx + J u^^^ipda^ 


0+5—21 \7„, |2—a 


|«+M 0 


|Vu| dx F CWuW^aF Cs 1 “ ll'^llLa([/) 


(4.10) 


FCs ^'^\\u\\1t'^j^^F j 


Define 


T(f) = ^ ^ f H{\'Vu{x,t)\)dx F [ u^^^FdcrF f u°‘dx. (4-11) 

2 Ju Jr Ju 


Adding p3.4p to p4.10p yields 
A + 1 


A 


[ ^{u^)‘fdx F ^£{t) F (Cl - 2e) [ 

Ju dt Jjj 

< C(1 + ||( 


|o+5-2|vv, ,|2-a 


\Vuy-^dx 


L°°(r) J l + ll L°‘{U) “T II “II L°‘{U) 


+ C\\u\\'^l ^ 

_l_ C{s~^^ F II 


+ C(e i-“ + 

2 —a 


2 — 0 , 
«— II 1 — 0. 


lL°°(r)^ll^llL“((7) 


IQ+MO 


llL°°(r) ^ll^llL“(!7) 

Choosing e sufficiently small such that Ci — 2e > 0, we derive 


|0+/i,l 


+ J l+tl (4.12) 


A + 1 
A 


f ^{u^)^dx F—£{t) < C{1 F 
Ju dt 


+ C(1 + 


2-a 
_— 11 1 —a 


ii,“(r)^iCiiL“({/) 


|O+M0 


lL°°(r))ll'U,|lL“([/) + ^11'“ 

2-0 
iL^(r) 

2 —a ^ 2—a 


a+S—2 

L^(U) 


+ C(1 + 


- II (l-a)(l-e) UL, ||“+Ml I 


J 


^-1 da. 


Note that a + (5 — 2<a<a + ;Uo<« + +i and 1 < We apply Young’s 

inequality for each norm on the right-hand side and obtain 


A + 1 
A 


f ^{u^)‘fdxF—£(t)<C f |+t| + C(1 + 

Ju dt Jy 


-|| (l-a)(l-e) 

llL”(r) 




Let t G (0, T). Integrating both sides of previous inequality in t, we obtain 


A + 1 


r / 

Jo Ju 


"^dxdr F 


A + 1 


H{\Vu{x,t)\)dx F f u^'dx 
Ju 


< 


£{0) - f u^~^^ipdaFC f f + C f (1 + ||+ llioc^r)^ )(1 + 

t/r t/Q'^r t/Q 


For the first integral on the right-hand side, applying inequality (|2.6I) in Lemma l 2 .1 1 with a = s = 
A + 1 = 2 — (5 and p = 2 — a, we have 


L 


— I u^~^^ipda < 


u^+^da 


< 


.|e j |Vu|2-“dx + C j |u|^+^dx + Ce"T^ J 
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for any e > 0. Now, using H{\Vu\) > CdVnp “ — 1) from (|4.3I) and applying Young’s inequality 
to the last two integrals with a > A + 1 + /tq we obtain 


- J u^~^^ipda < Css\\(p llLoo(r) J {H{\Vu{x,t)\) + l)(ki 


+ 




Selecting e = 


we obtain 


4C8(1+||¥’ IU°°) 


A + 1 


f f ^{u^)^dxdT + 
Jo Ju 


A + 1 


lu 


H{\Vu{x,t)\)dx + - / u^'dx 
2 Ju 


<£{H)+C + C\\ip + CKl + ||(/9 ||Loo(r))i “||(/? ||L°°(r)}“ ^ ^ ''o 

lo + (1 + ll+"liro;j)""®0(l + \\u\\lt^u))dT 

which gives 

f [ v}~^{u^)‘ldxdT + f H{\Vu{x,t)\)dx 

Jo Ju Ju 


'U 

< 


c[£{^)+lCi{t) + j\l 


+ 11'^ IIlooJ)' '"’)(! +(4.13) 


where 


(2 — a)a 




+ 11+ (i)llL°o^ ' + / / \(pt\<^-^-^dadT. 
Jo Jr 


Note that £1(0) < C£q and, by Young’s inequality, JCi{t) < CIC{t). Hence, we obtain from (|4.13|) 

f [ u^~^{u^)1dxdT + f H{\'Vu{x,t)\)dx 

Jo Ju Ju 

< C£o + CJC(t) + cl\l + + lklli+S)))dT. 

Last, using relation between H{x,t) and |Vu| again in ()4.14D . we obtain ()4.5I1 . 


(4.14) 


□ 


Now, we combine Theorem 13.31 and Proposition 14.11 to obtain a gradient estimate in terms of 
initial and boundary data. 

Theorem 4.2. Assume \4.4^ - IfT>0 satisfies 1^3.14 ) then for all t E (0,r], one has 

[ [ u^-^\{u%\‘^dxdT+ [ \Vu{x,t)\^-^dx<C£o + C}C{t), (4.15) 

Jo Ju Ju 

where £q and IC{t) are defined in Proposition \4-l\ and C > 0 depends on a, 

Proof. Using (j3.15jl to estimate the L“-norm of u in (j4.14p . we get 

rt 


[ [ u^-^{u^)UxdT+ [ |Vn|2-“dx < CTo + C/C(i) 

Jo Ju Ju 


a+n pt 


2-a 
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By using ()3.14p to bound the last integral, we obtain 





0 Ju 


{u^)^dxdT + 


lu 


\Vu\‘^-'^dx <C£o + CJC{t) + c(l + j u^ix)dx 


Then (I4.15p follows. 


□ 


5 Interior estimates 

In this section, we estimate the L°°-norm of the solution in the interior of the domain by using 
Moser’s iteration. 

Lemma 5.1. Assume a>2 — 5. Suppose Br and Bp, with R > p > 0, are two concentric balls in 
a compact subset of U. If T > T 2 > Ti > 0 then 


sup / u°‘{x,t)dx + 

ts[T2,T] J Bp 


< CaS, 


'T 2 JBp 


where 


S = 


u°‘dxdt + 


/Ti JBr 


u'^dxdt 


Tl JBr 


= C„(p, ii, Tl, Ts, T) C9«2 (i + I^^IT) (1 + 
with cg > 1 independent ofU, a, p, R, Ti, T 2 , T. 


a + S-2 


+ 


T 2 - Tl {R- p) 


2-a / ’ 


(5.1) 

(5.2) 

(5.3) 


Proof. Let C{x,t) = ^i(|a;|)C 2 (i) be the the cut-off function which is 1 on Tp x (T 2 ,T) and has 
compact support in Br x (Ti,T). More specifically, ^i(|x|),^ 2 (i) £ [0,1] and satisfy 


ei(N) = 


1 if jxj < p, 


0 if jxj > R, 
Also, there is C > 0 such that 


and f 2 {t) = 


0 if 0 < t < Tl, 
1 if T 2 < t < T. 


16! < 7 ^^ and |V^|< ^ 


T 2 -T 1 


R — p 


(5.4) 


In the calculations within this proof, notation C denotes a generic constant independent of a, p, 
Ri Tl, T 2 , T. 

Recall that <5 = 1 — A. Multiplying the PDE in (I1.26P by integrating over U, and 

using integration by parts, we obtain 


A [ u^-^f^dx= [ V ■ {K{\Vu\)Vu)u°+^-^fdx 


lu 


dt 


lu 


-(a-A) [ A:(|Vm|)|Vu| 

Ju 


2^a+5-2c2 


Cdx - 2 


[ K{\Vu\)V 

Ju 




Using properties (|3.2p . resp., (|3.1I) of function K{-) in the first, resp., second integral on the right- 
hand side of the last identity, we find 




A / u^-^i^^dx<-d^{a-\) I \Vu\'^-‘^u^+^-^fdx + d3 {a-X) I u^+^-'^tdx 


n+5-2c2^ 


lu 


dt 


lu 


lu 


+ 2d2 [ |Vu|^-“|u|“+^-ie|VC|dx. 

Ju 


(5.5) 
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Let e > 0. Applying Young’s inequality to the last integral of (15.5p . for conjugate exponents 
and 2 — a, we have 


a—1 I |^|a;+(5—a(^a|x-7tr|2— 

fu 


A f < —d^{a — \) f |Vm|^ + ^ 3(0 — A) f 

Ju ot Ju Ju 

+ e [ |Vu| 2 -“|M|"+^- 2 ^ 2 ^x + Ce“-^ [ 

Ju Jv 

Choosing s = then have 

L 

+ C{a-\) [ u»+^-^edx+ I 

Ju (« - Aji^ “1 Ju 

Then integrating the inequality in time from 0 to t gives 

d3(a - A) / f |y^|2-a^a+5-2^2^^^^ 

Jo Ju 


^a+5-2^2dx 


C\^^\^-^dx. 


A^ 

a dt 


u^i^dx - — [ u^^^tdx < ^ 

a Itt 2 Itj 


C 


q+( 5—a Aa I v-7 A12—a 




« Ju 


u°‘{x,t)^^{x,t)dx + 


< 


— [ [ u°‘^\Ct\dxdT + C{a - X) [ [ 

<a Jo Ju Jo Ju 


u^+^-^fdxdT 


+ 


C 


(a - A)(i-“) Jo 


\u 


a.-\-5—ac(i\x-7/:\2—a 


(5.6) 


u 


Using (15.411 to estimate and on the right-hand side of (j5.6p . the bound and support of ^{x, t), 
we have 


J — sup f u°‘{x,t)dx+ 
“ te[T2,T] J Bp 


d^ia - A) 


IT 2 JBp 




< 


AC 


f [ dxdr + C {a — \) 

J Tt J Br 


a{T2 - Ti) Jrp^ Jbj^ 


Note that A < 1 and a — A > 1. We then have 




'Ti JBr 

C 


u^+^-'^dxdr 

rT r 


{a-xy °-{R- py “ Jt, Jbr 


U^+^-^dxdT. 


J < 


c 


f [ dxdr + C{a — X) f f ‘^dxdr + 

J Ti J Bji J Ti J Bji 


c 


T 2 -T 1 

Inequality ()2.4[) gives + u°‘, hence 


{R-pf- Jt, 


\a-\-S—a 


dxdr. 


'Bp 


J <C{a- X) 


1 + 


1 




T 2 - Ti {R- p) 


2—a 


Ji, 


where Ji = u°‘dxdT + ‘^dxdr. Therefore, 


sup / u°^{x,t)dx < Ca^ 
t&[T2,T]JBp 


. 1 1 
1 + 7^-^ + 


T 2 - Ti {R- p) 


2—a 


Ju 


(5.7) 
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'T2 JBp 


\Vu\‘^-^u^+^-^dxdT < C 


1 1 
1 + ;^-^ + 


Ta-Ti {R-p)^-\ 


■Ji. 


(5.8) 


Applying Holder’s inequality to the second integral of Ji using exponents and we have 

Ji<{l + {\BR\Tf-^)S < 2(1 + \Br\T)S. (5.9) 

Hence estimate (15.11) follows from (|5.7I) . (|5.8I1 and (|5.9I1 . □ 

Proposition 5.2. Assume Ii2.29\) . Let Br, Bp, and T, T 2 , Ti and Ca be as in Lemma \5.1[ Then 

(5.10) 


(5.11) 

(5.12) 


where k is defined in \2.31\} . 

/ N def 1 x 0 ^ 

r = ryot) = a + o — 2, s = 5 ( 0 ) = 


a + 6 — a 

Aa = Cio(l + /3"^)^"“a®"“(l + IBrIT)^ (^1 + + ^J^_p^ 2 -a) ’ 

with cio > 1 independent ofU, a, p, R, Ti, T 2 , T. 

Proof. Applying Sobolev inequality (|2.32p to the ball Bp in place of Br, we have 

1-0', i 


J 


|l® • sup ( I | 

' ^ V / r? / ) 


te[T2,T] ''Jb, 


IT2 JBp 

where c = C 6 (l + exponent 9 is defined in (12.311) . and 

rT 


(5.13) 


I = 


\u 


a+S—a^ 


T 2 JBp 


dxdt + f [ ^|Vrt|^ “dxdt') 

Jt 2 JBp ’ 


Using (EU), we have 


/ < ^4 


T 


T 2 JBp 


I afi- 5 —a 


dxdt\ +CJ [ [ 

^ JT 2 'J Bp 


ocfi-6—a 


where C 4 = 2“+^-“. Then applying Holder’s inequality to the first integral on the right-hand side 
yields 

I <C4{T\Bp\)^^ [ [ lul^^dxdt + cJ [ [ |u|"+'^-2|y^|2-a^^^A 
JT2 JBp ^ ''T- In / 


T2 JBp 


Estimating the second integral on the right-hand side of (15.1411 by ()5.ip , and combining with (I5.13P 
give 


J < c^[{c^{T\Bp\)mS + C4(C„5)^) (C„5)1-®}“ 

11/ a — 5 a \ 0-|-l — 04 

-|(^C4(r|Hp|)^^+^5 + C4(C„5)^J+^+C„5j 1“ 


< c 

< 


c^{(C4(r|Hp|)^ +Ca)5 + C4C„“+^-“5re|' 
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Note by definition (15.211 that S = where ya = y \u\^dxdtj Thus, we find 

that 


J < c« 


{(C4(r|i?p|)^ +c„)(y" + 2/S+''"') + C4C“+^-“(y2 + 2/2+'-')^}' 

.,2 


(q: + 5 —2 )a 


< 


c^{(Miy“ + Miy^+^-^ + } 


where 




Ml = C4(T|i?p|) “+^-“ + Ca, M2 = C^Ca — <^400- 

—, then by ( 12 .4p we have 


Note that a + 5-2< < « < ^ 


(Q+(5_2)a 

a a+S-a ^ a-\-S-2 

UoLi UoL ^ iJa 


+ yr"““. 


Therefore, we obtain 


J< [3c«(Mi+M2)]^(y; + y^)^. 

Because < 1, a+'^s-a < 2, Cq, > 1 + |Bk|T > 1, and 1 < 6*4 < 4, we have 

Ml + M 2 < C*4(l + |i?/j|T) + Cq + C'IC'a — + Ca + C'fC'a ^ 3C*|C^ < 3(4 Cq)^. 

Also, c« < c. Combining (j5.16p and (j5.17p . we obtain 

J< [122cC2]^(y; + y^)i 


(5.15) 


(5.16) 

(5.17) 


Then (|5.10p follows. 


□ 


Iterating relation ()5.10p . we obtain the following local estimate for u. 

Theorem 5.3. Assume oq > 0 such that a = ao satisfies 112. 2 9\) . Let Br, with R > 0, be a ball in 
a compact subset of U, and T > 0, u G (0,1). Then 


u 


dBii/ 2 X{crT,T)) < Ci?,T,cr max I I 


0 (Br X (o,r)) > 11 “ 11 L“o (s^ X (o,r)) 


}. 


(5.18) 


where 


T = ll 

j=0 


— 2 + (5 

j ’ 
aoKi 


V = 


aoK% 


n , , 

j=o + 0 — a 


(5.19) 


CR,T,a — 2^^Cioao “(1 + -R ^)^(1 + |ilii;|T)^^l + + j^2-a) 


(5.20) 


with K* = K{ao) defined in 112. 3 1\) . constant cio as in Provosition I5.M and some positive number lv 
depending on oq. 


Proof. For j = 0,1, 2,..., let 

tj — (^^{1 “ Pi = ^(1 + ^)) Qj — ^Pj ^ 

where Bp. is the ball of radius pj having the same center as B^. Note 

(tT ./? 

tj - tj+i = ^: Pj - Pj+1 = ^> ^.1™ tj = , ^lim Pj = R/2. 
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Also, K* > 1. Let Oj = aoni. Then aj > oq gives K,{aj) > K{ao) = k*. 

Define Yj = ||w||l“j(( 3 ^.)- Note that (Qj)^o ^ sequence of nested cylinders. By Holder’s 
inequality we obtain 

Hence applying (IS.lOp to a = aj, p = Pj+i, R = pj, T 2 = tj+i and Ti = tj gives 


1 1 

j j 


< (i + |Qo|)^‘“^'A„; 

Using definitions in (I5.11|) and (I5.12|) . we denote 

Vj = r(aj) = Uj — 2 + 6 = oqkI — 2 + 5, 

>1 

Aj = (1 + IQol) ■ 


Sj = s{aj) = c?J{aj + 5 — a) = / {olqk,{ + 5 — a), 


Then 


y,+i < a;^' [Yp + Yp) . 


J V J 

We estimate Aj. We have from (I5.12h that 

A,< 


(5.21) 

(5.22) 


^ 2-^+^ / 2 .?+ 2 \ 2 -a', 2 

cio(l + IQol) • (1 + 2 /ii) a- “(l + IQol) } 

< 4ci„{l + R-‘)=(l + |Q„|)=(a„K;)'-“lfP«(l + T + tA-)" < A’l+, 


aT R^-^' 


where 


A 


R^T,a = max 


|i6k® “jd^CioOo “(1 + -R ^)^(1 + |<3ol)^(^l + ^ + j^2-a) }• 

1 n2 


Since k* G ( 1 , 2 ), we actually have 

A_r,t,o-= d^cioag °‘{1 + R ^)^(1 + |Qo|)^(^l + ^ + ^ 2 -a) — ^ 


(5.23) 


Therefore, 


i+i 

13 + 1 + ^R,T,(7 ^ ^ 


Vj > 0. 


(5.24) 


Since k* > 1, we clearly have + l)/'aj converges to a positive number. Note also that 


0 < In — = ln(l + 
• ^ ^ i ■ ^ 


a — 5 


< 


E 


a — 6 


j =0 


j =0 


aoni + 6 — a j^q aoni + 5 — a 


< 00, 


and 


0 < — In ^ In — = ln(l + 


2 + 5 




2 + 5 


aoni — 2 — 5 ao^i — 2 — 5 


< 00. 


(5.25) 


(5.26) 


j =0 ■' j=0 j=0 

Therefore, IF+f^^rj/oj) and n^Q(sj/aj) converge to positive numbers p and u, resp., given by 

(CTI) . 

By (15.241) . and applying Lemma [A.21 to sequence (1^ )^q, we obtain 


\u 


\L^{Ba/^x{aT,T)) = 3 ^ {‘^^R,T,aT ^ayi{Y^,Y q], 


for some positive number U3. Then the desired estimate (15.181) follows. 


□ 
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Remark 5.4. Inequality i5.18\) obviously leads to the quasi-homogeneous estimate il.24\ ), which 
was proved in for equation 

ut = V ■ A{x, t, It, Vu) (5-27) 

with the homogeneous structure 

A(x, i, tt, Vtt) • Vw > c|nP“^|Vii|^, \A{x,t,u,Vu)\ < (5.28) 

Due to the non-homogeneity of function K{-)> see our equation \1.21\) cannot he converted 

to {5.21^ , i5.28\) . Therefore, above inequality 115.18\} is an extension of ^L2^ to the class of 
equations lil.21\) with non-homogeneous structure LI.1\) . 

Now, we bound the L°°-norm of u, in any compact subsets of U, in terms of the initial and 
boundary data. 

Theorem 5.5. Let U' he an open, relatively compact subset ofU, and a = ao satisfy i2.29\} . Then 
for T > 0, and 0 < e < min{l,T}, one has 

ll'^lli/°°([/'x(£,r)) ^ (1 + T) j"’ (5.29) 

where oj,p.,u are the same as in Theorem \5.,‘A 

In particular, ifT>0 satisfies I13.12\) for a = ao, then 


u 


\L°°{U'x{e,T)) < Ce ^‘^(l + r)^“ max 


\3iJ . 


u 


ao 




V^,ao 


CXQ 


}, (5.30) 


where 


U, 


ao 


(t) = {{^ + 


h‘'l,aQ 

ao 


-C: 


3,00 


pt 2— 

(1 + ||^-(T)||ioo(r))(T^^)^dr}"^. (5.31) 

(5.32) 


Furthermore, ifT>0 satisfies ^3.14^ for a = uq then 

\\u\\L^{U'x{e,T)) < -I-r)3‘^+!^/“0(']^ _|_ ||tto H^ao ([/) )'^- 

Above, constant C > Q depends on U, U' and oq- 

Proof. Denote R = ^dist{U',dU) > 0. Because the set U' is compact, there exists finitely many 
Xi E U', i = 1,..., m for some m E N such that is an open covering of U'. For each 

i, we have from (|5.18l) with e = ctT that 


W\\L’^{Bji/2(xi)x{£,T)) — CR,T,e™ax ||ri||^Q,Q^(yx(o,T))|! (5.33) 


where 


CR,T,e = [2}^cioal-%l + \BR\Tf{l + R-^)\l + e"! + i?“-2)2}" < CRe-^^{l + Tf^ 

for some positive number Cr depending on R and oq. Summing up the estimates ()5.33p in i, we 
obtain (|5.29l) . 

In case T > 0 satisfies (I3.12p for a = oq, we use (I3.13P to estimate the L“°-norm in (|5.29l) . and 
obtain 


\l°°{U'x(£,t)) ^ Ce 2“(l + T^‘^) max I f / f \u\^°dxdt\°‘° , (^ f f \u\°'°dxdt^ 

hQ J Bpi JO J Br 

T , Ji_ 


0 J Bji 


< Ce- 2‘^(1 + T^) max { ( / ““ }, 


This proves (15.301) . 

In case T satisfies (j3.14p for a = oQ) we use estimate (I3.15|) in (I5.29p and note that v > p,. Then 
we obtain (15.320 . □ 
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6 Global L'^-estimates 


In this section, we estimate the L°°-norm on U of the solution u(x, t). We perform Moser’s iteration 
on the entire domain, and take into account the effect of the Robin boundary condition. Thanks 
to the contribution of the boundary terms, calculations need to be much more meticulous. 

Lemma 6.1. Assume a > max{2 — <5, n/io}- Let 6, pi, fi 2 , be defined as in LemmalKM 

//T > T 2 > Ti > 0 then 


sup / n“(ix + j 
t&[T2,T]Ju Jto Jv 


T 2 JU 


with 


S = 


iTi JU 


1 2—a^^a+S—2 

T 

Ti JU 




‘^dxdt < MS 

( 6 . 1 ) 

. a+S-2 

^^dxdtj , 

( 6 . 2 ) 

-|- E 3 + E 4 + E 5 ), 

(6.3) 


where constant cn > 0 is independent of a, Ti, T 2 , T, while 


fj-l —S+2 


^i = (T 2 -ri)-'(|R|r)'^+Mi, E 2 = {\U\T) , i?3 = (|C/|T)“+«||(^-|Uoc(rx(o,T)), 


Ml(a+Mn) Ml-MO 2 ii£ fi. 


with d^ > 0 defined in IRUI) . and 


def 


F4: = FA,a = F2 + i = 


2 — a 


{l-a){l-er 


(6.4) 


Proof. Let ^(t) be a smooth cut-off function with f{t) E [0,1], ^(t) = 0 for 0 < t < Ti, (^(t) = 1 for 
T 2 < t < T, and 

161 <ci/(r2-ri), (6.5) 

for some C > 0 independent of Ti,T 2 . 

Multiplying the PDE in (jl.26j) by integrating over U, and using integration by parts 

yield 

A J u^-^f^dx = J ^ J V-(iL(|Vu|)Vu)u“+^-i,e^ 


= -(a-A) [ K{\Vu\)\Vu\\^+^-^fdx- [ ^u^fda. 

Ju Jr 

Using relation ()3.2p we obtain 

A [ u^-^i^^dx <-dfia-X) [ \Vu\‘^-V+^-^fdx + d3 {a-X) [ + J(t), 

Ju ot Ju Ju 


where 


Since a — A > 1, we have 


J{t) = \W {t)\\L^{v) j^iL^{x,t)f{t)da. 


A /" u" ^C^^dx < —ds [ |Vu|^ ‘^f^'^dx + dsa [ ^dx + J(t), 

Ju ot Ju Ju 
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Using the product rule on the left-hand side of the inequality we have 


—-f- [ u^^'^dx + d^ [ |Vm|^ ^^‘^dx<— [ u°‘^^tdx + dsa [ dx + J{t). (6.6) 

a dt Ju Jjj a 


Integrating from 0 to t, and taking supremum for t G [0, T] give 


A 

“ [o,r] 


sup f u^{x,t)^\t)dx + d 3 f [ < — f f u^C\Ct\dxdt 

[o,ri Ju Jo Ju “ Jo Ju 

11 

Jo Ju 


+ 2d3a 

Using the trace inequality (|2.25|1 and noting that ^{t) is independent of x, we can estimate J 


^a+S 2 f J(^t)dt. (6.7) 

40 


by 


J{t) < \\g) (t)||2,oo('r)|2e j '^C'^dx + J u°'^'^dx + 

a g+Hi 

u'^fdx'^ “ }, 


-D3,„e i-“ 


a a + ^O 

u^^'^dxj “ 


where fj ,2 is defined by (j2.26p . Hence, 


1 f J{t)dt 

Jo 


< 


W llL°°(rx(0,T))“ 
+ 1 - 


\4e [ [ |u|“+^-2|Vup-“^2^xdT + (7 [ [ u^i^dxdr 

^ Jo Ju Jo Ju 

/ T r q:+ho pT r 

(^J “ dr + 2U»4,ae"'^' J [j u^^fdx'^ “ dr}. 


/o ^JU 

Next, applying Holder’s inequality to the last three integrals yields 


2 

where 


/ J{t)dt < Ae\\ip ||L°°(rx(o,T)) 
Jo 


\u 


q:+( 5—21 V7,i, 12—a 


^0 Ju 


|Vu| ^idxdT + C\\ip ||L<x=(rx(o,r))>/o, (6.8) 


a 1 ^l(a+Mn) n-W 01 + p.Q 

Jo = (|t/|r) “+'"1 y “+'‘1 + “(“+Mi) r a+,ii Y°‘+n ^ D 4 ^a£~'^^\U\~Y 


with Y = Jjj u°^^^^^’^dxdT. Combining (16.71) and (16.81) with properties of ^{t) gives 
-sup/ u“(x,t)^^(t)dr + (da - 4e||v7"||ioo(^x(0,T))) / / |Vup"“u“+'^"^C^dxdr 

“ [o,r] Ju Jo Ju 


< 


C 


< 


a(T2 — T] 
C 


vj:i 


u°‘^dxdT -|- Idsa 


I 


u 


0+5-2a2 


0 JU 


^ dxdr + C\\ip ||Loo(rx(o,T))do 


—^ + 2 q: + (5 —2 


a{T2 - T] 


_(|[/|T)-+MiFa+Mi +2d3a(|t/|r) “+^^1 Y +C\\ip \\l^(tx(o,t))Jo- (6.9) 


Choosing e = 


_^3_ 

8||r“llLoo(ax(0,T)) 


, and using properties of ^(t), we have 


- sup [ u“(x,t)dx + ^ [ [ |Vu|2-“u“+^-2dxdr 

“ te[r 2 ,T] Ju 2 Jt 2 Ju 

ot q; + 5 —2 

< C{Ei/a + H3)y “+M1 CaE2Y “+''i + CE^Y^^+i^i + CE^Y. 
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" HpF T' 

Note from the choice of the cut-off function ^{t) that Y <Y = Then 

sup f u°‘(x,t)dx < CaJi and f f \Vu\‘^~°‘u°‘^^~‘^dxdT < CJi, 
te[T2,T]Ju Jt2 Ju 

— at ^ ol-\-5 —2 Q~l~M0 

where Ji = {Ei/a + E^)Y + aE 2 Y “+'"i + £^4^“+''! + E^Y. 

Comparing the powers of Y in Ji’s formula and applying (|2.4I) , we have 

Ji < 3{Ei/a + aE 2 + £3 + £4 + E^)<S < 3a(£i + £2 -|- £3 + £4 + E^)S. 


( 6 . 10 ) 


Hence, 


sup f u°'{x,t)dx+ f [ |Vttp “^dxdrYCaJi, 

t&[T2,T] Ju Jt2 Ju 

and we obtain ()6.ip . 

Proposition 6.2. Assume a > max{2 — 5, repo}- If T > T 2 > Ti > 0 then 


\W\\l'^°‘{Ux{T2,T)) < (^lklll,a+Mi(;7x(Ti,r)) + II^IIl“+'‘i (C/x (ri,r)) j ’ 

where k is defined by 112. 31\) . pi is defined by 112. 26\} . 

s def a{a + 5- 2) , . def a{a + pi) 

r = r{a) = -;-, s = s{a) = 


a + gi 


a + 5 — 


(M1+a-i5)a a g 

Aa = cuot + A4“+''i + 

with A4 defined by 116. 31) . and C12 > 1 independent of a, Ti, T 2 , T. 

Proof. We use Sobolev inequality (I2.30p in Lemma [231 

rT 




where c = the numbers 9 and k are defined in ()2.3ip and 


I = 


u 


"+'^-“dxdt+ / / \u\^+^-^\Vu\'^-^dxdt 


IT 2 JU 

Applying inequality (I2.3p , we find that 

rT 


IT2 JU 


q+ 5 —a 


I<CA 


( f / +<^5( J |u|“+^-2|Vup-“dx(it) (6.15) 


□ 

( 6 . 11 ) 

( 6 . 12 ) 

(6.13) 


O pT p 1 C \ 0 ^ 

/ \u\'^°'dxdt\'^°‘ <c+Tiil^ ■ sup ( / \u\^dx') (6-14) 

T2 Ju J t&[T2,T]^Ju J ’ 


where £5 = 2 “+<5-“ ^ = 2“+^-“ . Applying Holder’s inequality to the first integral on the right-hand 
side of (j6.15p with conjugate exponents and , we get 


I < 


(Ml+a-<5)a / fT f 

C^{T\U\)++n)++^-<^)yj J |u|“+^idxdtj 


1+/J1 


T 
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Next, we use (16.11) to estimate right-hand side of (I6.16p . Hence combining ()6.14p and (16.1611 yields 




1 r r tMl+n-OJa , a _ q, 1 U - i n'l ' 

J<c-{ C'5(r|[/|)(“+'^i)(“+^-“)5“+w -hC'5(A45)^^+^ 

1 r r (Ml+n-<^)a _ a ~ a ~~\ 

<c^| C'5(r|t/|) (“+'‘1 )(“+*-“) 5“+''! -I-C'5(A45)“+^-“ -|-Cq,5 

1 r (Ml+a-i)a; ~ _a_ _ g „ g 'I ^ 

= S C*5(T|C/|) (“ + 4l)(a + «-a) ^a + iii a + S-a a + S-a K _ 

~ ~ Q; ~ q: 

Since < 1 < , we use (|2.4I) to estimate 5 < 5“+''i -|- 5Thus, we have 

1 r (Ml+a-'^)a a , a a - g -| ^ 

J < c«a [6*5 (r| [/| ) (a + fil)(a + 5-a) _|_ a+Ml J^a+Ml -|- [1 -|- C^]M a + — a a + 5 — a s 


Denote Ji = {Jx^ Jjj\u\°‘~^^^dxdt)°‘+n . Then by definition p6.2l) of S, and applying (12. 2 p . resp. 
(12.31) . we find 

a(a+5 — 2 ) 

“+"1 , 


resp., 


Therefore, 


where 




0!{oi-\-fj,l) (q: + 5 —2)o 


- a(Q' + <5 —2) a(a+fii) (a: + 5 —2)a .. A 

J < c-{MiJf + Ml “+'^1 + M2 Ji"^^+^ + M 2 J 1 “+'■“ }“ , 


(Ml+a-'5)a g _ g 

Ml = C'5(r|C/|)(“+'^iH“+^-“) +A4a+Mi^ M 2 = C5(1-fC5)A4^^+^ 


(6.17) 


Since f < a, < s, we apply (12.41) to estimate the first and last summands on the right-hand 

side of (16.171) by 

(q: + 5 —2)cx Q:(q:+^ —2) Q;(Q: + /ij^) 

ja J oi-\-5 — a ^ J a+Ml _|_ J a+5 —a 


Then it follows 


J < c^|3(Mi + M2)(Jr 



(6.18) 


Since a > 2 — 5 > 2{a — 5), we have C 5 < 2 , hence Mi -|- M 2 < 97Vi. Note also that c« < c, then 
we obtain (16.110 from (I6.18p . □ 


Now to perform the iteration we need to start with an initial exponent K(ao)ao such that 
K(ao)ao > ao + d'i,ao: that is, K(ao)ao/(«o + A*i,oo) > define 

( 6 . 

The following properties are useful in later iterations. 


K{a) = 


K{a)a 

01 + /ri,a 


Lemma 6.3. For a € ((2 — a)a!*/(l — a), 00 ), the functions a —>■ ^i,a in ^2.26j> and a —>■ ^ 4,0 in 
are decreasing, while the function a hi{a) in 116.19\) is increasing. 
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Proof. First of all, we note that K{a) defined in (|2.3ip is increasing, while 9a defined in (12.241) is 
decreasing. Since fii defined in terms of 9 in (I2.26p is increasing in 9, and 9 = 9a is decreasing in 
a, then, as a composition, pi^a is decreasing in a. Similar argument applies to g.A,a- 
Next, if a' > a we have K{a') > K{a) and pi^a' < 9 - 1 ,a, hence 

K(a')a' K{a)a' K{a)a' K(a)a 

K{a j = — -> —-> —-> -= K[a). 

OC + 9x^0.' OC + 9\^a' OC + Pi^a Ot + pi^a 

Therefore, K{a) is increasing in a. □ 


We construct two sequences of exponents in order to implement Moser’s iteration. (Regarding 
the notation, the numbers a^’s below are newly constructed and are not the exponents in (jl.lOl) .) 

Lemma 6.4 (Construction of a^-’s and /3j’s.). Let 

2 + 02-a)(2 + i)-l 


Assume oq > (l + x*)a*, let 9^ = 9ao, 9* = 9iao, n* = ^(ao), and k* = ^(ao). Define the sequence 

iPj)T=o by 

fio = ao + //*, fij = Rifio for j > 1. (6.21) 

Then: 


(i) K* > 1 . 

(ii) There exists a strictly increasing sequence such that 

= “j + 91,aj Vj > 0. (6.22) 


(hi) For all j > 0, 

9aj < 9^, 9i^aj < 9*, L{aj) > R^. (6.23) 

(iv) For all j > 0, 

aj < L-iPo, (6.24) 

and there exists a number k* > 1 such that 


Oij > h{ao Vj > 0. 


(6.25) 


Proof, (i) Using definitions of k, pi and po in (I2.3ip . (j2.26p and (j2.23p . the inequality k* > 1 is 
rewritten explicitly as 


1 1 

ao( 1 + (a - (^)(- 

a* ao 


fr| + Gaoid - 


> ao + 


1 - 


«o 


which is equivalent to 


a* (1 - a)(l - 6»ao)’ 


Using formula (I2.24h for 9ao , we convert this inequality to a quadratic inequality in uq as 


(1 — a)aQ — 2(2 — a)a*ao + (2 — a)a* > 0. 
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Its positive solutions are 


, 1 + y/2 — a. 

ao > (IH-^-)«*• 

1 — a 


(6.26) 


This is satisfied by our choice of ao- Therefore /t* > 1. 

(ii) It follows from definition of (dj that {/3j)^Q, is unique and strictly increasing. Consider the 
equation 

(6.27) 


def 


I3j =x + = f{x) 


where, for x > 0, 




I — 0(x) ' (1 — a)(I — 0(x)) ’ 

with 0(x) = We have f{x) = 1 - ■ Then f{x) > 0 if 

l/0(x) > 1 + that is, X > a*|+ ij. 

y/n 11 — o V y/n / J 

Note that < y/2 — a, we already have from (I6.26|) that 

ao > 


1 ,2-a 

-(— 1=^ + 1) + 1 

1 — a y/n 


a*. 


Hence, / is strictly increasing on [ao,oo), /(ao) = /3o and /(oo) = oo. Since the sequence (/3j) is 
strictly increasing, we have for any j > 1 that jdj > /3o, and hence the number aj = f~^ildj) solves 
(I6.22p . Clearly, the sequence (aj)^Q is also strictly increasing. 

(iii) By the monotonicity of 9a, k{a) (Lemma 16.31 and its proof), and the fact 

Oij > ao, we have 

daj<0ao=9^, n{aj) > K{ao) = K{aj) > R{ao) = R,, 

(iv) From (|6.22p . we have aj < fdj and hence inequality (I6.24h follows. 

By (16.2211 and the fact gLi^a is decreasing in a, we have for j > I that 


_ _ l^j ^ T /^l,ao ^ I d'O 

\ - “T 


1 + 


oo/q;» —1 


I3j-1 


aj-i 


Oij—i ao 1 — 


Set ao = (1 + x)a*, then x > x*. Thus, 


K* < 


Hence, 


We have 


(l-a)(ao/“*-l) 

2 — a 

aj-i ' a*[(l —a)x —1] ckj-i ' n[(I — a)x — 1] 

2 — a 


a,- 


+ 


//o(l - a) 


a, 


+ 


a 


7 .. def _ 

— l>£ = fC* — 1 — 


Cf 7 —1 


n[(l — a)x — I] 


(6.28) 


K*ao _ ^ _ (1 + (g - 5){l/a^ - I/ao))ao _ ^ _ 


ao + /i* 


* 


(a — 5)x — /r* 

(1 + x) + /i* 


K* — 1 = 


ao + Rt.: 










































Generalized Forchheimer Bows of isentropic gases 


29 


Note that hence, 

^ ho(l + ^ao(!-«)) ^ (g - ^)(1 + ^ (g - 5){x + 1) 

1-Oao (l-g)x-l ■ 


Thus, 


We aim at finding x such that 


K* — 1 = 


X 

1 + X 


2—a 


+ 


1 

(l—a)x—l 

1 

(l—a)x—l 


e = 


X 

1+x 


— + 
2-a ^ 


1 

( 1 — a )£—1 
1 

(1—a)x—1 


2 — g 

n[(l — g)x — 1] 


> 0. 


(6.29) 


If inequality ()6.29p holds true then we choose k* = 1+e > 1, and by (j6.28l) . ajjaj-i > 1+e = k*. 
Thus, aj > K*Q!j_i, and by induction, aj > itiao for all j > 0. 

It remains to verify (|6.29p . For x > inequality (16.291) is equivalent to 


2 — g,^ „ 2 — g 


(1 — ayS:'^ — 4(1 — a)x^ + (1 + 2g- )x + 2-> 0. 


n 


n 


(6.30) 


Since 2 — > 0, a sufficient condition for (|6.30l) is 


(1 — afx^ — 4(1 — a)x + (1 + 2g — —^) > 0. 


2 — g. 


n 


Solving this inequality gives 


X > 


2 + ^/3-2g + ^ 
1 — g 


= X* 


which is satisfied by the choice of x. The proof is complete. 


□ 


The final preparation for Moser’s iteration is to estimate in ()6.1ip . 

Lemma 6.5. Let ao be a positive number such that a = ao satisfies 112. 29\) . Then one has for any 
OL ckQ that 

Aa < Ca^^{l + r)^®(l + _ 2^^)^(1 fo Wl’ lli“(rx(o,T)))^^) (6.31) 

where > 0 and C > 0 depend on ao but not on a. 

Proof. In this proof, C denotes a generic positive constant depending on ao, but not on a. 

In order to estimate A^, we estimate M. from (16.31) first. Note that 0 < //2 < F2,ao^ hence 
(d 3 / 4)“^2 < (fi Then from (16.3p . we find 

M < Ca^il + ^r^)(l + \\T-\\L^irxio,T))r^{i\U\T)^i + (If/IT)"^ 

Ml(a+Mn) IJ-l-m Ml I 

+ D3,a\U\ “(“+"1) T “+'"1 +D4^a\U\~\. 
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Above, we used the fact that > 1 is the maximum among the exponents of \\ip ||Loo(rx(o,r))- 
Next, using definitions of and in (\2.27\\ and (12.281) . we have 


1 


Hi Ml —^+2 

a+Mi + {\U\T) “+'"1 


M < Ca^{l + ^)(1 + \\ig> ||L°°(rx(o,r)))^'‘ 

a/ , r N (2-a)(l + g(l-a)) 2-a (l-a)(c«+fio)S Ml (a+Mfl) Ml-MO 

+ 2®("+^-“)c* a~\U\ -^— -Pl “(“+'' 1 ) 


e{a + S-a) (2-a)(l + e(l-a)) ^1 'I 

+ 2 1 -e pa) \Up 

where constants are as in Lemma For convenience in calculations below, we keep using c* to 
denote max{l, c*}. Also, using 9 < 1 and m < a, we have 


M < Ca^il + 


^ I e(a + S-a) (2-a)(l + e(l-a)) 

, 2/1 , )2^^Pa) (!-)(!-«) (l + ||<^-||Loc(rx( 0 ,T)))'‘" 


T 2 - Ti ■ 

Ml 


r Ml mi-' 5+2 (l-a)(a+Mo)^ , MlCo^ + Mo) Ml-MO ui '» 

■ U\U\T)^+n + {\U\T) <-+n + \U\ “ +-(-+mi) T “+w + |[/|” [. 

Set 0* = 9ao and /x* = /ri.ao- Since a > ao, we have da <9*- Hence 

(2 - a )2 


( 2 -q)( 1 + 6 >( 1 -q)) ( 2 -a)^ ( 2 -a)^ 

(c*a) < (c*a) < (c*a) where zi = 


(l-a)(l-0*)’ 


6{oL-\-5—a) _ 

Next, we want to bound 2 by some number independent of a. From (I2.24D . 


0(a + (f — a) = 


a + 6 — a 


(1 — a)(a/a* — 1 ) ’ 

which, due to the fact that a* > a — <5, is decreasing in a. Hence, 

Hpf 6(oi+5—a) Z 2 

9{a + 6 — a) < Z 2 = 0*(ao -|- <5 — a) and 2 i-® < 2i-®» 


For exponents of |H|r, we note that 


Hi Hi ~ ^ d- 2 ^ ^ 
a + fii' a + /ri “ 

For the remaining power of T, 

Hi - Ho ^ ^ 

a + /ii “ 

For the remaining powers of |f7|, 

/I + /^o ^ ^ r, Hip + Ho) ^ Hi ^ 

{l-a)e - <- <2, —-—^< H*, 

a a a[a + a 

Also, we have for the power of ||v^~||Lo°(rx(o,T))- 

def 2 Cl 

Hi — 2^3 Hi,ao T^j r7:j ) 

due to the decrease of in a, see Lemma [631 
So we find 

M < + ||92-||Loo(rx(o,T)))"^(l + \Up+^pl + T). 
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Summing up, we obtain 


M<M = + ll¥^"llL-(rx(o,T)))"^ 


Ta-Ti 

Since A4 > 1 and < 2, using (j6.32l) in (|6.13p we have 

n (mi +a-S)a _ 

Aa < ci 2 a^"“[(r|[/|)(“+M)(“+' 5 -“) + 2M^]. 

Note that a — 5 = (1 — a)/io and q: + 5 — a>2 — a, hence 

(iui + a - 6)a _ /io(2 - a)a ^ /io(2 -a) ^ Fo 


< 


(6.32) 


(6.33) 


(6.34) 


(6.35) 


(a +/ii)(a + (5 - a) (1 - 0)(a + pi)(a + S - a) (1 - 0*)(a + 5 - a) 1-0* 

Hence by (j6.32l) . (16.331) and (I6.34p . 

ia < C'a2-“+2(2+^i)(i + _^)2(1 + + ||<^"||Loo(rx(0,r)))""^- 

Hence we obtain p6.3ip from p6.35p with = 6 — a + 2zi, pQ = max {2, j^^}, and ^7 = 22 : 3 . □ 

Applying iteration process, we obtain: 

Theorem 6.6. Assume 

ao > max{2 — 5, (1 + x*)a*} (6.36) 

with X* defined by i6.20\) . There are C, jl,D,uJi,uj 2 ,uj 3 > 0 such that ifT>0 and a G (0,1) then 

1 , 


\mL^{Ux{crT,T)) ^ 


where /3o = oo + hi,ao ■ 


<C 1 + 


aTJ 


(1 + r)^2(l + ||(/7 ||L°o(rx(0,T))) 


UJ3 


• |ll'“llL/3o(f/x(0,T))’ ll'“llL'3o(f/x(o,r)) 


}, (6.37) 


Proof. Note that a* >a — 6 and npo = a*(2 — a)/(l — a), then it is easy to check that oq > 2(a — 5) 
and ao > npo. 

Let /is, /ig and /i 7 be defined as in Lemma [631 and aj, fij, k*, it*, it* be as in Lemma 
For j > 0, let tj = cjT( 1 — i), Qj = U x {tj,T), and define Yj = .)■ 

Applying (|6.1ip of Theorem 16.21 with a = Oj, T 2 = tj+i and Ti = tj, we have 




:r{aj) 


s(aj) 


T)) 


= AZi Ik 


'LG{Uxitj,T)) 


+ k 


LPi(Ux{tj,T)) 


where fj = f (aj) and Sj = s{aj), see formula (16.1211 . 

By part (hi) of Lemma [6~il K{aj)aj = K{aj)/3j > = fij+i, then by Holder’s inequality 


Combining the above two inequalities give 


Yj+i < Aj^ {Yp + Yp)‘^i , where Aj = |Q^+i| 


(6.38) 
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Now we estimate Aj. From (|6.31l) . (|6.38l) . the fact that aj < /3j+i, and (16.241) . we have 

2 i+i 


Aj < C(1 + (1 + +rr(i + 

1 


lL°°(rx(0,T))) 


U7 


< C(1 + |Qo|)(«e/3or4^(l + ^)2(1 + rr(i + 


where 


= max{4«f ,C/3o^«(l + —f(l + |[/|T )(1 + rr(i + 


lL°°(rx(0,T)))^^} > 1- 


Hence 


i+1 


From (I6.25P we have 




OO . ^ 1 • I 1 

V^<-V^<oo. 

Q/ . 


(6.39) 


i=i 


“0 — Ki 


Note that 


1 > — = 

a, 


Q^j ~t“ (5 — 2 Q^j “t“ (5 — 2 k^^oiq “t“ (5 — 2 


> 


> 


4 " CXj + d-* 


> 1 - 


1 < ^ = 


+ Ml,a 


■j < “i + M* 


< 


Kf ao + M* 
K^ao + M* 


= 1 + 


M* 4“ 2 — 6 

'•j ’ 

K^ao 

M* 4“ ® ^ 


Oi a,■ + 5 — a Oj + 6 — a kiao + 5 — a kiao + S — a 


X.J '-V 

Then it is elementary, see (I5.25h and (15.261) . to show that the products 


r, _ A 


o-j T — 2 


aj + Mi,aj 


and k = n°T, 


aj 4- Mi.Oj 
aj + 6 — a 


converge to positive numbers. By (I6.39D and Lemma lA.21 we obtain 

limsupF,- < ( 2 Hr,(T,<^)'^ max{l)f, 
j^oo 


(6.40) 

(6.41) 

(6.42) 

(6.43) 


where uj = G YlT=i ^ with G = YlT=ii^k/ak) E (0, oo). 

Note that 

(2^t,o-,vp)‘^ < ^*(1 + 4- T)‘^^(l + \\ip ||L°°(rx(o,T)))^ 

where oji = 2uj, W 2 = (1 4- Me)*^ and ws = jiju. Then estimate (|6.37l) follows (I6.43F 




□ 


Remark 6.7. (i) The exponents jl and D in ( (g.37D are given by ^6.4^ but can, in fact, be replaced 
by simpler and more explicit ones such as 


- _ TTOo - 2 + 6 

M - ■ 


M = n-0- 


aoki + M* 


aoki + M* " aoki + 6 — a 

Indeed, it follows from estimates and Ijd-fl]) that fi < ft < k < k, and then applying \2.4^ 

gives 

max{y(f, To"} < < 2 maxjy^f, Tq"}- 

(ii) Estimate \6.3'll\) is a global version of the improvement il.24\ ) on interior estimates. See 
also m for a similar global result for degenerate equations with the use of De Giorgi’s iteration 
instead. 
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We now have global L°°-estimates in terms of the initial and boundary data. 


Theorem 6.8. Assume oq satisfies 116.36\) . Let /3o = oq + Then there are positive numbers 

C,uji,uj 2 ,uj 3 , such that: 

(i) //T > 0 satisfies i 3.1 Si) with a = fio, and 0 < e < min{l,T}, then 

\\u\\L°°{Ux{e,T)) ^ Cs {1 + {I + \\ip ||Loo(rx(0,r)))‘^® 


• max 


{(I Upfit)dty\(^j Upfit)dt)y, (6.44) 


\l°°{Tx{o,t)))‘^^ ■ (6.45) 


where Ufifit) is defined by I15.31\) with fio replacing uq. 

(ii) IfT > 0 satisfies Ii3.14\ ) with a = fio, then 

\W\\l°°{Ux{6,t)) < Ce ‘^1(1 + (l + ||mo|Il/3o(;7)) (1 + II 

Proof, (i) Applying (|6.37l] to aT = e > 0, we have 

ll^llL°°({/x(£,r)) < C{I + Py^{I + ||¥?“||L°°(rx(o,T)))‘^® 

yj j \u\^° dxdy° , j J \u\^° dxdty°. (6.46) 


max • 


Using (j3.13p we have 


lu 


\u{x,t)ydx < Ugfit). 


(6.47) 


Then (I6.44p follows (j6.46p and (j6.47p . 

(ii) If T > 0 satisfies (|3.14p with a = fio, then (I3.15P gives 


J u^°{x,t)dx < 2 ^ + J UQ°{x)d3i 


Combining this with (16.461) yields 

ll'*^llL°°(i7x(£,r)) < + Ty^{l + ||<^“||Loo(rx(o,r)))‘^® 


• max 


J (1 + J Uq° {x)dx)dty°, (2 y (1 + J u^°{x)dx)dty°'^. 


Since H > fl, then we obtain (I6.45|) . 


□ 


A Appendix 

Let {yfi^Q be a sequence of non-negative numbers. 

Lemma A.l. Let (aj)^g and be sequences of positive numbers with 

00 00 

d = < 00 and fi = JJ fij exists and belongs to E (0,oo). 

j=o j=o 
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Suppose there is A>1 such that 


yj>0. 


(A.l) 


Then {yj)^Q is a bounded sequence. More specifically, for all j > 1, 


y, < -ft-.. 


(A.2) 


and consequently. 


limsupyj < A^^y^, 

j^oo 


(A.3) 


where Bj = max{l, l3rnfim+ifim+2 ■ ■ ■ fin M < m < n < j}, and B = limsup^^go Bj. 

Proof. Applying (|A.ip recursively, we have 

y,+i < 

<^ * * * <^ ^oij+aj_i/3jH-Qj_2/3j_i/3j+Q:j_3/3j_2/3j_i/3j+...H-Qo/3i/32.../3j ^ 

It follows that yj+i < Hence, we obtain (jA.2jl . Taking the limit superior of 

(jA.2p as j —>■ oo, we obtain pA.3p . Note that H < oo by Cauchy’s criterion. □ 

The next lemma is the main adaptation used in this paper. 

Lemma A.2. Let Kj > 0, Sj > > 0 and u)j > 1 for all j > 0. Suppose there is A > 1 such that 


1 

— r,- , Si 


yj +1 < A "2 {yd + y.^) -j Vj > 0. 
Denote fij = rj/nj and jj = Sj/nj. Assume 


CO 

d = — < oo and the products 7 j converge to positive numbers fi,'y, resp. 

i=o j=o 


OO OO 


- Kj 

j=0 ^ 


Thcfi 

yj < ( 2 A)‘^^“max|yJ°-'^^"\yJ°-^^"^| Vj 
where Gj = max{l, 7 m 7 m+i .. .^n 'H < m < n < j}. Consequently, 


> 1 , 


(A.4) 


limsupi/j < ( 2 A)^“ max{yg, i/q }, where G = limsupGj. 


,7 


(A.5) 


J^OO 


J^OO 


7' 1 

— r,' , Si 


Proof. We prove (jA.4p first. Define a sequence {zfi^Q by zq = yo and Zj+i = {zfi + zfi)'"^ for 
j > 0. Then 

yj < H J — O' ("^'6) 

Therefore, it suffices to bound zj. We consider three cases. 

Case 1 : zq > 1. Clearly, Zj > 1 for all j. Hence, 


Zj+i < A^j (2z-0 < (2A) 7 z! 


7i 
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Then using Lemma lA.ll we have 




(A.7) 


Together with (IA. 6 j) . we obtain (jA.4p . 
Case 2 : zj < 1 for all j > 0. Then 


iti J_ Li fl 

Zj+i < A'^j {2z]^)'^i < {2A) "i z^ 


Applying Lemma lA.ll gives 




(A. 8 ) 


where Bj = max{l, I3mf3m+i ■ ■ ■ l^n ■ 1 < n < j} < oo. Note that Bj < Gj. Then, again, (lA.Sp 

and ()A. 6 P yield (|A.4p . 

Case 3 : There is jo > 1 such that zj < 1 for 0 < j < jo and Zj^ > 1. Applying (IA. 8 I) to 

1 < j < jo, 

Zj < {2A)^i < ^2Afi Y.\=l>^il^i _ (A.9) 


Same as Case 1, Zj > 1 for all j > jo- Then applying (IA.7D for j > jo gives 




Jo 


(A.IO) 


where Gj^j = supjl, 7jo+m7jo +m+i ■ ■■'Ijo+n ■ f < rn < n < j - jo} < oo. Using inequality (IA.9P 
with j = jo to estimate zj^ in (lA.lOp , we have for j > jo that 




Since zo < 1, A < 'Ji for all i, and Gjp 7 jp 7 ^ 0 + 1 ... 7 j_i < Gj, we obtain 

Zj < (2A)™'^^'f'^'’'0’J^TOO'joO'jo+i-7j-i}Ei7o< ^2A)^^ ^i=o (ATI) 

for all j > jo- Then, (|A.9p . (jA.lip and ()A. 6 p imply (|A.4I) . This completes the proof of ()A.4ll for all 
cases. 

Now, taking the limit superior of ()A.4I1 as j ^ 00 yields ()A.5D . □ 

Note that the numbers G,d,j3 and 7 in inequality (jA.5h are explicitly defined. 
Acknowledgment. L. H. acknowledges the support by NSF grant DMS-1412796. 
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